Bose— Hubbard model for universal quantum walk-based computation 
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We present a novel scheme for universal quantum computation based on spinless interacting 
bosonic quantum walkers on a piecewise-constant graph, described by the two-dimensional Bose- 
Hubbard model. Arbitrary X and Z rotations are constructed, as well as an entangling two-qubit 
CPHASE gate and a SWAP gate. Quantum information is encoded in the positions of the walkers on 
the graph, as in previous quantum walk-based proposals for universal quantum computation, though 
in contrast to prior schemes this proposal requires a number of vertices only linear in the number of 
encoded qubits. It allows single-qubit measurements to be performed in a straightforward manner 
with localized operators, and can make use of existing quantum error correcting codes either directly 
within the universal gate set provided, or by extending the lattice to a third dimension. We present 
an intuitive example of a logical encoding to implement the seven-qubit Steane code. Finally, an 
implementation in terms of ultracold atoms in optical lattices is suggested. 



I. INTRODUCTION 

Quantum walks have proven to be a fruitful alterna- 
tive to the quantum circuit model for the construction 
and description of quantum information processing tasks 
[l| . The framework they provide has allowed for the con- 
struction of efficient quantum algorithms 0, including 
some previously unknown under other models l3|-|8| and 
alternative formulations of Grover's search M as well 
as reproductions of other known results (Tol. Ill|. Fur- 
thermore, it has since been shown that quantum walks 
are universal for quantum computation. Three distinct 
models have been described, based on continuous-time 
quantum walks [l2j, discrete-time quantum walks [l3| . 
and what we refer to as the discontinuous walk, a hy- 
brid method in which a continuous-time quantum walker 
takes discrete steps across a set of time- varying graphs, 
undergoing perfect state transfer through a subset of the 
graph at each stage [l3|. In each case, the use of a single 
walker to encode n qubits leads to a graph with f2(2") 
vertices. 

One straightforward adaptation that allows for the ex- 
ponentially growing Hilbert space required to encode n 
qubits without necessitating a similarly sized graph is to 
employ multiple quantum walkers. As a step in that di- 
rection, a framework for the description of multiple non- 
interacting discrete-time quantum walkers, distinguish- 
able or not, has been put forth [l^. It is straightforward 
to see though that if there are multiple walkers on a graph 
yet they do not interact, then there is no meaningful 
difference from the situation of multiple distinct copies 
of the single-walker situation. Non-interacting walkers 
evolve independently of each other, so while they can be 
useful if many runs are required to build up statistics of 
the output state, no new dynamics can be present. 

In the discrete-time quantum walk, effective interac- 
tions between otherwise non-interacting walkers can be 
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induced through the sharing or swapping of coins jl6| . 
and entanglement can arise between a single walker and 
its coin [i7|. Continuous-time quantum walks on the 
other hand have no coin degree of freedom, so must rely 
on inter-walker interactions to introduce additional dy- 
namics and generate entanglement. There is evidence 
that on a given graph, two interacting continuous-time 
walkers are more computationally powerful than either a 
single walker or two non-interactin g o nes when applied 
to the graph isomorphism problem [l8| . 

For the case of the continuous-time quantum walk, 
the Bose-Hubbard model [l^ provides a natural method 
with which to describe interacting bosonic walkers, and 
its relationship to a single continuous-time quantum 
walker has been discussed in the context of particle trans- 
port and entanglement generation [20| . The model is 
well studied in the realm of condensed-matter physics; 
in particular it describes strongly correlated systems well 
[21| , providing an excellent description of features such as 
the transition between the superfluid and Mott-insulator 
phases, as well as havin g be en implemented and well con- 
trolled experimentally |22l |. Furthermore the standard 
Bose-Hubbard Hamiltonian does not address internal de- 
grees of freedom, allowing us to explicitly consider the 
encoding of quantum information within position states, 
in the spirit of quantum walks, rather than a more con- 
ventional encoding in spin states. 

Standard quantum walk-based schemes for quantum 
computation employ a position-based encoding. This 
results in a spatially delocalized qubit under circum- 
stances where the quantum walk formalism can be di- 
rectly mapped to the circuit model. In this work, we 
assign a pair of vertices to each bosonic walker, corre- 
sponding directly to the computational states |0) and |1). 
Gates are implemented by making instantaneous changes 
to the graph supporting the walkers, and we explicitly re- 
quire couplings to doubly occupied states in which two 
walkers interact on one vertex, despite the fact that such 
states do not encode qubits. The graph is fixed except at 
the instants of change, so continuous dynamical control 
is not required. The main idea is to design a discrete 
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sequence of graphs, such that at the beginning and end, 
the spatial separation of the walkers is maintained with 
one walker per pair of vertices, yet during their evolution 
the walkers have the opportunity to interact whenever 
two walkers occupy one vertex. This allows states with 
doubly occupied vertices to be harnessed as an advan- 
tage in constructing logical gates, rather than having to 
be strictly avoided as a source of decoherence during gate 
operations as in previous proposals f23j . 

In this work we provide an important correspondence 
between the discontinuous evolution of interacting in- 
distinguishable particles on polynomial-sized graphs and 
that of a single discontinuous quantum walker on an ex- 
ponentially larger graph. That is, a number of walkers 
linear in the number of qubits to be simulated n, walk- 
ing on a vertex set of a size also linear in n, has the 
same power as a quantum walk with a constant number 
of walkers (in particular, one walker) on a graph with a 
number of vertices exponential in n. 

The structure of this article is as follows. In Sec. |TT] 
we briefly review the Bose-Hubbard Hamiltonian, then 
define an encoding of computational basis states within 
a subset of the position states available to the bosons 
in the system it governs, and describe a set of physi- 
cal operations that result in a universal set of one- and 
two-qubit gates on the computational space. We then 
describe in Sec. IIIII how single-qubit measurements can 
be made on encoded qubits, and provide an example 
scheme for adding a layer of quantum error correction. 
Finally in Sec. II VI we discuss the possibility of a proof-of- 
concept implementation of our scheme, before providing 
some concluding remarks in Sec. |Vl 



II. POSITIONS AS COMPUTATIONAL STATES 
A. The Bose— Hubbard Hamiltonian 

We use the Bose-Hubbard model to describe multiple 
interacting continuous-time quantum walkers on a graph 
as spinless bosons hopping on a lattice, with on-site in- 
teractions. The graph G is defined by a set of vertices 
V and associated edges E (- V x V . The vertices corre- 
spond to the lattice sites, while the edges indicate allowed 
tunnelings. The Bose-Hubbard Hamiltonian is 

?i = -J^cIc^.-f-^^n,(n,-l). (1) 

Here c| creates a particle on vertex i and Ui = cjcj is 
the number operator on said vertex. The first sum runs 
over neighboring vertices such that £ is an 
edge of the lattice graph. J is the tunneling amplitude 
between adjacent sites, and U is the on-site interaction 
strength between the particles. 

We refer to the graph G, on which multiple quantum 
walkers appear, as the primary graph. This primary 




i=l i=2 1=3 



FIG. 1. Cartoon representation of six vertices confining n = 3 
particles whose positions encode a computational state of the 
form \'4>) = |0) ® |1) ® The physical Hilbert space H for 
this setup is 56-dimensional. The Hamiltonians we prescribe 
couple the encoded 8-dimensional computational space C to 
a total of 24 of the 56 basis states. 

graph and the n walkers evolving on it comprise a phys- 
ical system that can be described by a set of states and 
allowed transitions between pairs of them. One can con- 
sider each of these states as corresponding with a vertex 
in a secondary graph, in which each allowed transition 
appears as an edge. The evolution of n walkers on the 
primary graph then maps exactly onto the evolution of 
a single walker on the secondary graph. Note that the 
number of vertices in the secondary graph is at least ex- 
ponentially larger than the number of vertices in G. 

We make use of n bosonic quantum walkers to encode 
n qubits. On a primary graph of 2n vertices we assign a 
pair of vertices to each qubit, and for simplicity visualize 
the vertices as being arranged in a 2 x n grid structure; 
see Fig. [T] for a cartoon example with n = 3. Label the 
vertices Vi^x S V , with i S {1, . . . , n} and x G {0, 1}, and 
let \Qi_x) be the vacuum state on vertex Vi^x in row x of 
column i. The vacuum state of the system is then 

|vac) = |Oi,oOi4---0„,oO„,i), (2) 

which when convenient can also be expressed as 

|vac) = |00)i---|00)„, (3) 

where the first (second) entry in each pair is understood 
to correspond with row (1) of the array. A walker is 
created on vertex Vi^x by the operator c|^^. 

The physical Hilbert space H of n bosons confined to 
such a lattice is larger than the 2"-dimensional compu- 
tational Hilbert space C which their positions encode. 
The encoding is accomplished with a subset He C H, 
such that I He I = 2", along with an isomorphism enc : 
He — )■ C. Generally we can use C and He interchange- 
ably without ambiguity. In order to generate entangled 
states in C we make use of a Hamiltonian that couples He 
to states outside of the computational space, i.e. those 
in Hj^ = H \ He. Under our scheme, a physical state 
1 4*) e He provides a valid encoding of a computational 
state on n qubits, |4')e = enc(|4')), while any state 
|a) G H that has non-zero support in Hj^ has no such 
encoding. Since the computationally entangling Hamil- 
tonian couples He to H_l, the evolution of the system 
restricted to C appears to be non-unitary in that the 
magnitude of the projection of the physical state of the 
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system onto the subspace He can be different from unity. 
However, we show that an arbitrary initial state with 
support entirely in He will, at well-defined times in its 
evolution under the Hamiltonians we specify, map onto a 
state that also has no support in Hj^. That is, all proba- 
bility will return to C and the physical system will once 
again encode a valid computational state. 

The ability to construct an entangling Hamiltonian un- 
der this encoding relies on the indistinguishability of the 
bosons involved as the encoding we use requires the op- 
erations of interchanging two particles on the lattice and 
of swapping them twice, returning them to their initial 
configuration, to be identical. This is crucial for preserv- 
ing the mapping from the physical system to the com- 
putational space that requires one bosonic walker to be 
localized to two vertices of the primary graph. 

We assume that unless an explicit operation is being 
performed, the lattice is deep enough to prevent all tun- 
neling. In this default configuration the system Hamilto- 
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This results in a primary graph that is completely dis- 
connected — 2n vertices with no edges, E ~ 0. Each 
two-site column of the lattice represents a qubit, with 
the upper and lower sites corresponding to the |0) and 
\1) computational basis states of that qubit, respectively. 
An n-particlc state | represents a valid computational 
state if and only if 



(5) 



for each i from 1 to n, i.e. if there is one particle in each 
column. Physical states that satisfy this criterion are 
mapped onto computational ones in a canonical way. 



|0), ^ c]JOO) 



and 11). ocJilOO) 



|10)., 
101),. 



(6a) 
(6b) 



We use underlines to indicate encoded computational 
states, in contrast with physical Fock-number states. The 
system is initialized in the state 
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On 
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This is a zero-energy eigenstatc of "Ho, and as such is sta- 
tionary while the lattice is maintained and the primary 
graph has no edges. 

The dimension of H is (^^) , in general much larger than 
|C| = 2". This may seem inefficient or even detrimental 
at first sight, but it is this access to a larger Hilbert 
space during evolution that allows for the implementa- 
tion of a two-qubit controlled gate through only on-site 



interactions with no internal degrees of freedom. We im- 
pose Hamiltonians that couple He to Hj^ in such a way 
that at well-defined times the resulting unitary operators 
block diagonalize as C/ = C/c ffi U±. That is, all popula- 
tion initially in a valid computational state returns to 
the computational subspace at the end of the evolution, 
despite having been transferred through the larger space 
at intermediate times. 

Before discussing our two-qubit gates, in particular an 
entangling CPHASE gate and a nearest-neighbor swap, we 
first provide a simple method for implementing arbitrary 
single-qubit operations. 



B. Single-qubit operations 

Consider the ith encoded qubit, with computational 
basis states |0)i and \V)i encoded according to ([S]) in a 
physical state j^") satisfying the single-particle condition 
dni). The single-qubit operations we construct below triv- 
ially preserve this condition, as they do not couple \'^) 
to states outside of the computational space. 

We first show how to implement an arbitrary X ro- 
tation on qubit i. By lowering the height of the lat- 
tice barrier between the two sites we can set the system 
Hamiltonian to 



n 



X,i 



-J 



H.c. 
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where Jx,i depends on the height of the lowered barrier 
between the sites. The primary graph in this case is 
composed of 2(n — 1) disconnected vertices and one copy 
of the connected two- vertex graph, on the vertices 
encoding qubit i. The action of 'Kx,i on the basis states 
([6|) is simply 



nx,^■■\\G)J ^ -(5,,Jx,,|01)j, 

|01), ^ -%Jx,z|10)„ 



(9a) 
(9b) 



where (5.^ is the Kronecker delta. In the computational 
space this acts as an X operator on qubit i. 



|1), ^ -JxM^ 



(10) 

(11) 



and does nothing to the other qubits encoded in |\1/). 
The unitary operator generated by evolution under this 
Hamiltonian for a time t is then Ux,i(i) = Rx.ii— '2Jx.it), 
an X rotation of the ith qubit. 

Next we construct a Z rotation on qubit i. To do so, 
we maintain the initial height of the barrier between the 
two sites while applying a local potential to the site. 
The physical Hamiltonian becomes 



n 



z.inis 



Ho 



(12) 



and corresponds to a primary graph of 2n disconnected 
vertices, with a self-loop attached to the second of the 
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FIG. 2. An n- walker graph on 2?! vertices, on which appro- 
priately initiahzed bosonic walkers will undergo X and Z ro- 
tations on the encoded qubits i and j, respectively. 



if the applied-potential-to-hopping ratio is tuned to be 
Vh^i/Jh,! = 2. Not only does this one-step process re- 
quire fewer dynamical controls, but for equal hopping 
and interactions terms (i.e. taking JH,i = Jx,i and 
Vz,i = Vha) results in a run-time that is an order of 
magnitude shorter. This example is unlikely to be the 
only such shortcut to additional gates available by judi- 
cious choices of further Hamiltonians. 



two vertices that encode qubit i. The Hamiltonian acts 
only on the computational basis state, as 



nz^^ ■■ 11), ^ -S^JVzAl)J, 



(13) 



SO the resulting unitary is Uz.iit) = e"^^^-^*^^ Rz.iiVz,it) 
— a, Z rotation of qubit i, up to an unimportant overall 
phase. 

Given an angle 6 £ [0, 2tt) wc can therefore perform 
Rx{()) on qubit i by evolving under Hx.i for a time 



An -9 
'2Jx,i 



or e ' Rz{0) by evolving under T-Lz,! for a time 

tzM = 

Vz.i 



(14a) 



(14b) 



Given sufficient freedom in the ability to set the tunnel- 
ing rates and on-site potentials, it is possible to enact 
either of these single-qubit gates on each qubit simulta- 
neously, with different values of 9 on each one, and have 
them finish at the same time. That is, given a set of 
angles 9i and a choice of gates Oi G {RxiRz,!}, the 
values of Jx.i and Vz,i can be chosen such that in a fixed 
time t, the operations Oi{9i) are simultaneously applied 
across all qubits i S {1, . . . , n}. Fig. [5] shows an example 
multi-walker primary graph that applies an X rotation 
to qubit i and a Z rotation to qubit j . The combination 
of these operations allows for the execution of arbitrary 
single-qubit unitaries in three steps by decomposing the 
corresponding rotations of the Bloch sphere using Euler 
angles. 

Finally, we note that while these single-qubit opera- 
tors are sufficient to implement a Hadamard operation 
on qubit i in three steps, as 



m = Rx,^{^T|2)RzA^^|2)RxA^/2), 



(15) 



it is also possible to obtain a Hadamard in a single step 
with the Hamiltonian 



H.c. +Ho. (16) 



This simple approach of effectively turning on the Hamil- 
tonians for X and Z rotations simultaneously results in 
the application of a Hadamard gate on qubit i, up to an 
overall phase, at time 



tH,i — 



2\/2J}i.i 



(17) 



C. Generating entanglement 

We now show how to generate entanglement between 
adjacent qubits, in the form of a controlled phase gate, 
CPHASe(0). Wc must expand our discussion to a 2 x 2 
contiguous sub-block of the entire lattice, which we take 
to be the vertices vi^^ and for x S {0, 1}, and we 

again assume that the initial n-qubit state \^) satisfies 
condition (0), containing one bosonic walker per two-site 
column. Restricting our attention to two walkers on four 
sites, the physical space in question is 10 dimensional. 
Four basis states correspond to the computational basis, 
and we will couple these to an additional four physical 
states. The remaining two physical basis states can be 
ignored so long as the initial state is a computational 
one. For convenience we drop subscripts and write the 
Fock states on the four vertices in question as a single 
ket when there is no ambiguity from doing so, taking for 
example |0110) = |01)i|10)i+i. The two-qubit computa- 
tional space is encoded in the Fock number states 



4o4i.olOOOO) = 11010) ^ 

4o4i,i 10000) = 11001) ^ |01), 
4i4+i,ol0000) = 10110) ^ |10), 
4icl+i,il0000) = 10011) o 111). 



(18a) 
(18b) 
(18c) 
(18d) 



The six physical basis states |110G), |0011), |20G0), 
|0200), |0020), and |0002) have no computational inter- 
pretation. We can think about the two-qubit situation 
in two complementary ways. The first is as we have de- 
scribed, a four-vertex primary graph on which there are 
two interacting quantum walkers. The second is to treat 
each of the physical basis states as a vertex in a new 
single- walker secondary graph on 10 vertices, with edges 
prescribed by the transitions allowed under the Hamilto- 
nian in question. For example, the physical and compu- 
tational bases under the Hamiltonian we will introduce 
for generating entanglement arc represented in Fig. [31 
with the primary graph in (a) and the corresponding sec- 
ondary graph in (b). 

To implement a CPHASE gate, we decrease the barrier 
height between the sites corresponding to |l)i and |l)i-|-i 
so that the system Hamiltonian is 



He 



-Jr: 



H.c. 



(19) 



Let Ucp{t) = exp(— iHcpi) be the time-evolution operator 
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(a) 



(b) 



i i + 1 



|1010) 



|1001) loiio) 



Q 

2000) 

10020) 



loioi) 

10002) 



|1100) looii) 



10200) 



FIG. 3. (a) The primary two-walker subgraph on which en- 
coded qubits i and i + 1 undergo a CPHASE operation has four 
vertices and a single edge, (b) The corresponding secondary 
graph, describing the couplings among Fock states under the 
Hamiltonian Hqf, in columns i and i+1 of the physical graph, 
has 10 vertices to represent the 10-dimensional Fock space H 
of two bosons on four sites. Four of these are selected to 
encode the computational space He. Edges between vertices 
correspond to allowed transitions, and self loops to on-site 
interactions. 



generated by Hop, where we have dropped the H' sub- 
scripts when no ambiguity arises from doing so. Then 
for any time t, the |00) computational state evolves as 
;7cp(t)|1100) = 1 1100). The state encoding |01) couples 
outside of the computational basis, to |1100); at time t, 

t/cp(t)|1001) = cos(JcpO|1001) +«sin(Jcpt)|1100). (20) 

In order to guarantee that Ucp maps computational states 
to computational states, we must therefore require it to 
act for a time tcp,fc = kir/Jcp, < fc € Z. This also satis- 
fies the requirement that |10) be mapped to the compu- 
tational basis. Specifically, for lip) G {|1001), |0110)} we 

have UcAtcp,kM = 

The fourth computational state on two qubits, |11). 
evolves within a three-dimensional subspace spanned by 
10101), 10200), and |0002), as illustrated by the largest 
connected component of the secondary graph in Fig.[3jb). 
Determining the action of C/cp on |0101) shows that for 
a walker initially in this state to return entirely to the 
computational basis at time tcp,fc requires 



exp 



27riW16 



El 

'-"CP 



= 1. 



(21) 



This is accomplished if the ratio of the on-site repulsion 
to the hopping rate is tuned to be U/Jcp = ^/rn? — 16, 
for any integer m > 4. The final requirement to guar- 
antee that Ucp maps the computational space onto itself 
at time tcp.fc is that the product mk be even. We can 
guarantee this by setting k = 2, defining t^p = 2tt/ J^p. 
In this case the action of Ucp decomposes into its action 
on the computational basis states, and its action on the 
rest of Fock space. In the canonical ordering of the com- 
putational basis states, we have 

UcAtcp) = diag (l, 1, 1, e— ("+^/^^^)) e Ui, (22a) 

= CPHASE((pm) © f/J;, (22b) 




FIG. 4. Distribution of available phases satisfying (|23[) . 
Except for the two trivial values, <p4 and <p5, no phase is a 
rational multiple of vr. For large m the phase goes as —8n/m, 
modulo 27r, returning to as m — >■ oo, which corresponds to 
the limit of zero hopping. 



where tj^ acts only on the subspace orthogonal to the 
computational one, and 
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(23) 



This is a non-trivial entangling phase for any m > 5, and 
the available values are depicted graphically on the unit 
circle in Fig. [H 

If we instead choose an odd value for fc, including k = 1 
which results in a shorter run-time for the gate, the even 
values of m > 4 result in the same set of phases but 
indexed by to/ 2 rather than to. The resulting gate in 
this case is {Z ® Z)cPHASE((y3„/2). 



D. SWAP gate 

In order to provide a universal gate set with only a 
nearest-neighbor entangling gate, a two-qubit swap gate 
is also required. A variation of the X rotation performed 
by 'Hx,i allows for a straightforward implementation of 
this by the simultaneous lowering of the potential barrier 
between sites \x)i and |a;)i+i, for each x € {0,1}. The 
system Hamiltonian in this case is 



H.c, 



Ho. (24) 



This acts non-trivially on the four basis states satisfying 
(l5|), coupling the computational space to all six of the 
remaining physical basis states. As with the analysis of 
the CPHASE gate, we can place restrictions on the avail- 
able parameters such that the action of the operation 
restricted to the computational space is unitary at the 
end of the evolution. The action of T-Ls on |1001) and 
|0110) requires that we set the swap time to be 



ts,k 



{2k+l)TT 

2Js ' 



0< kez. 



(25a) 
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Under this restriction, the action of i-Ls on |1010) and 
|0101) further requires that 




2A; + 1< Z e Z. 



(25b) 



When these conditions are satisfied, the action of Us{t) = 
cyip{^iHst) at time t = ts.k block diagonaUzes such that 
its effect on the computational space is that of 



Usits,k) 



with 



fe-'°"' 

0-10 
0-10 

V e-'""" ' 



= 1+ y'P - 4k{k + l) -1. 



(26) 



(27) 



In general this provides us with a second entangling gate, 
unless a is an integer. In that case, if a is even then (|26p 
is equivalent to {Z (E) Z)swAP, and if it is odd then we 
obtain the gate —swap. 

The problem of finding values for k and / that satisfy 
the conditions p5|l while resulting in an integral value of 
a reduces to finding Pythagorean triples, of which there 
are infinitely many. One such combination, which re- 
sults in the minimal time for the operation, is fc = 1 and 
I ~ 5. This pair leads to the parameters tg = 3TT/{2Jg) 
and U/Js = 16/3, and implements a two-qubit nearest- 
neighbor SWAP gate, up to an overall phase, as required 
for a universal gate set. 



E. Connection to discontinuous walks 




FIG. 5. (a) Hopping parameters turned on to enact gates with 
two qubits present, (b) The corresponding graphs encoding 
couplings among the ten states of H are presented in gray 
for each operation, with the four computational states of He 
highlighted in light gray. 



(a) 



Rx (X> SWAP 



I 



(1) (2) (3) 




FIG. 6. When three qubits are encoded, (a) three edges 
among the six vertices of the primary graph implement an 
Rx <8> SWAP gate, (b) The corresponding secondary graph 
is much larger, and is given by the cartesian product of the 
one- and two-qubit graphs for Rx and SWAP. Note that those 
vertices which are not connected to the computational states 
have been omitted in an attempt at clarity. 



Implementing a sequence of computational gates re- 
quires that we toggle a set of potentials on and off in a 
prescribed order, affecting the evolution of a set of quan- 
tum walkers in the process. This is in the same spirit as 
the single-walker discontinuous walk we have previously 
proposed [l^ , though that work makes use of a number 
of vertices exponential in the number of encoded qubits 
whereas in the current multi-walker scheme only linear 
growth is required in the number of vertices. 

As in the prior models of universal computation by 
quantum walk [l2l . [isj . the single- walker discontinuous 
model employs a set of 'rails' for encoding the quan- 
tum information. These rails are linear graphs that al- 
low propagation of information from left to right, inter- 
spersed with so-called 'widgets' that enact transforma- 
tions on the state of the walker as it passes them. There 
is one rail for each computational basis state, which leads 
to the exponential growth in vertices. In Fig. [S] we show 
how a time sequence of two-walker primary graphs on 
four vertices corresponds to a single discontinuous walker 
on a larger secondary graph, and draw a connection to 
the rail model. Since 2n = 2" for ri = 2 the number 



of rails is equal to the number of vertices in this case, 
though in Fig. [6] we show what a single step of a discon- 
tinuous walk on three qubits looks like, on six vertices 
in the primary multi-walker case and on eight rails in a 
56-dimensional space for a single walker on the resulting 
secondary graph. 

III. MEASUREMENTS AND ERROR 
CORRECTION 

Under a single-walker model for computation, a mea- 
surement of the computational states of m < n qubits 
can be accomplished either in a single step with a set 
of fi(2"') physical measurement operators, each spatially 
delocalized over 57(2""™) vertices, or as a sequence of m 
single-qubit measurements, each requiring two measure- 
ment operators delocalized over 17(2""^) vertices. As- 
sume the best-case scenario of there being a one-to-one 
mapping between the computational basis states and a 
(sub)set of 2" vertices. For x e {0, . . . , 2"* - 1}, let 
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X = Xi ■ ■ ■ Xm be the m-bit binary expansion of x. To 
make a measurement of qubits {«fe}^i, let Vx be the set 
of integers between and 2" — 1 whose n-bit binary ex- 
pansions have bit value Xk at bit position i^. Then the 
measurement operators required are 

|p.= 5^|«)h| . (28) 

Clearly there are 2™ values of x, and for each x there 
are 2""™ unspecified bit values so \Vx\ = 2""™. To im- 
plement m' single-qubit measurements, simply set m = 1 
and repeat m' times. This exponential growth in the spa- 
tial extent of the required measurement operators is why 
previous quantum-walk models have been proposed pri- 
marily in terms of computational capability, and not in 
terms of possible physical implementations. It also makes 
the prospect of implementing quantum error correcting 
codes impractical at best. 

In contrast, the Bose-Hubbard-based multi- walker 
model presented here allows for the measurement of m 
qubits with localized measurement operators. Define 

= {cly\0^.^){Q^,^\{^S = \K^){b.J, (29) 

the projector onto the state with exactly b walkers on 
vertex Vi^x- The qubit encoded on vertices Vi^ and Vi^i 
can then be measured with the operators 

M,,o - P/.o' ® P}°i ® I^, (30a) 

M^^ = Pfo' ® P^^ ® Ii, (30b) 

where is the identity operator on every vertex not la- 
beled by i. These two operators form a measurement 
basis for the system since the allowed states when a mea- 
surement is to be performed have at most one walker on 
vertex Vi^x, so the identity operator on that vertex can 
be resolved as 

Ir.x - |0,,.)(0,,,| + = P^°^ + pI^. (31) 

Furthermore the probability to find either zero or two (or 
more) walkers on the pair of vertices {wj.o, Vi^i} vanishes, 
so all other two-site projectors 

« + (32) 

are not part of the basis for the states encoding qubit i. 

As in the single-walker case we can perform a single- 
qubit measurement with one pair of measurement oper- 
ators, but in this case each operator is spatially localized 
to two vertices regardless of n. Performing an m-qubit 
measurement in a single step, as opposed to by making 
m single-qubit measurements in succession, requires a set 
of 2™ measurement operators of the form (PU]) . Again, 
each one requires projectors only onto single vertices. 

With the ability to perform single-qubit measurements 
on the multi-walker system in a straightforward man- 
ner comes the ability to implement those quantum error 




FIG. 7. Addition of the seven-qubit Steane quantum error- 
correcting code to the Bose-Hubbard-based multi-walker 
scheme for universal quantum computation. The computation 
scheme described in Sec. Ull takes place on the horizontal row 
of light-gray vertices, in which the two in region A have been 
singled out as the qubit to be encoded in this cartoon. The 
14 vertices of region B (including those in ^4) provide seven 
physical qubits with which to encode a single logical qubit. 
Region C provides six additional ancillary qubits that can be 
entangled with those of the logical qubit in order to perform 
syndrome measurements as a sequence of single-qubit mea- 
surements. Note that the ordering of the |0) and |1) rows has 
been reversed in C with respect to B\ this allows a straight- 
forward method of generating entanglement between the two 
regions by way of the CPHASE gate that results from an edge 
between the |1) vertices of neighboring qubits. 

correcting codes (QECCs) that rely on encoding a sin- 
gle logical qubit in multiple physical one s l24| - |26j . One 
such code is the seven-qubit Steane code j27|. a so-called 
CSS code capable of detecting and correcting arbitrary 
single-qubit errors and of being implemented fault tol- 
erantly since operations on encoded logical qubits can 
be implemented by way of local operations on the un- 
derlying physical qubits. This and other QECCs can be 
implemented directly within the multi- walker framework 
already discussed in previous sections, constructing the 
requisite operators to within a desired error tolerance 
from a polynomial number of gates from the universal 
set provided. An alternative, which we now discuss for 
the seven-qubit Steane code in particular, is to use a 
straightforward extension of the qubit-encoding scheme. 
Consider extending the 2 x n grid into the third dimen- 
sion, to a 2 X 7 X n array of vertices or, when it is eas- 
ier to perform syndrome measurements as single-qubit 
measurements rather than the four-qubit measurement 
operators described by the code, to 4 x 7 x n, as de- 
picted in Fig. [T] Each pair of vertices encoding a single 
qubit in the original scheme described above, as in Re- 
gion A of Fig. [71 has a further six pairs extending from 
it along the third dimension of the lattice. The result- 
ing total of seven pairs at a given position i within the 
lattice. Region B, allows for the creation of a single log- 
ical qubit. Six further ancillary qubits. Region C, can 
be entangled with the physical qubits of the logical one 
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in order to implement syndrome measurements on single 
qubits. Making the vertices that correspond to the |1) 
states in the two regions adjacent allows entanglement 
to be generated between them by the CPHASE gate ([^^ 
with the straightforward addition of edges between the 
logical region and the ancillary one. 

Logical operations on these additional qubits can of 
course be performed in exactly the same manner as on 
the original qubits described in Sec. |IT1 Those gates are 
dependent on certain couplings between pairs of vertices, 
but are not inherently reliant on a grid- or lattice-like 
structure in the placement of those vertices; such a set- 
ting merely provides a useful visualization method, and 
yields an obvious tic to possible extant physical systems 
as we discuss in Scc. lIVI They can therefore all be applied 
between pairs of qubits in this new arrangement as well. 

Having seen that our multi-walker Bose-Hubbard- 
based approach to universal computation provides the 
advantage of being able to implement computational er- 
ror correction, previously unaddressed within quantum 
walk-based quantum computation, we also note that the 
dynamical nature of the underlying graph introduces a 
second potential source of error. If a gate is timed in- 
correctly, it is possible for the system to end up in a 
physical state that does not encode a computational one 
- besides being flipped and dephased, a qubit can be lost 
altogether. When this possibility is allowed for, the mea- 
surement operators (|30|) no longer sum to the identity 
and wc must introduce a third operator. 



(33) 



This provides a simple method for detecting the loss of 
a qubit, but not a mechanism for recovering from it. 
The possibility of losing the walker exists in the orig- 
inal proposal for universal quantum computation by a 
single quantum walker [l^ . and a similar timing issue 
is present in the discontinuous-walk scheme [14| as well 
as any potential physical scheme for the simulation of a 
discrete-time quantum walk with a continuously evolving 
physical system. These are not fundamental limitations 
on quantum walk-based computing, but must be consid- 
ered in any serious attempt to engineer such systems. 



IV. PHYSICAL IMPLEMENTATION 

A key feature required of a physical system that is 
to implement this proposal is the ability to provide the 
specified 2 x n lattice — or three-dimensional lattice, if 
that method of implementing a QECC is chosen — such 
that the tunneling amplitude between neighboring sites 
is close to zero when no gate is being enacted. This 
corresponds to the qubits' having a long coherence time. 
In practice this can be accomplished in either the 2D 
or 3D case by creating a three-dimensional lattice with 
isolated wells, and ignoring any vertices outside of the 
primary graph to be implemented. 



The additional requirements of state preparation, ma- 
nipulation, and read-out are common to any quantum 
computer [28| . and we discuss them here in the context 
of our proposal. State preparation is accomplished by 
the initialization of the jO)**" state, by loading one boson 
onto the first of each pair of sites in the primary graph. 
To manipulate the encoded qubits, it must be possible 
to selectively increase the tunneling amplitude between 
given adjacent sites in order to enact X, CPHASE, and 
SWAP gates, and to change on-site potentials to enact Z 
gates. Finally to read out the result of a computation it 
must be possible to measure the positions of the bosonic 
walkers within the lattice, as discussed in Sec. IIIII 

The experimental scheme proposed in Ref. [l^l includes 
many of the features required, though it makes use of 
adiabatic processes for gate executions. A combination 
of this approach with the sudden potential-landscape 
changes discussed in Ref. [l^ is more appropriate to the 
current scheme, and there have been significant experi- 
mental advances in the intervening years that offer the 
promise of a proof-of-principle implementation. 

There is an obvious connection between our scheme 
and optical-lattice experiments, and several experiments 
satisfy the above requirements. One option is to use a 
liquid-crystal display (LCD) as a spatial light modulator 
[30|. Such devices generate arrays of microtraps holo- 
graphically based on the pattern of opacity and trans- 
parency present on the easily programmed LCD screen. 
The traps have been used to store single neutral atoms 
per site, address individual sites, and measure the loca- 
tions of trapped atoms within the lattice [30|, thus pro- 
viding a means to create the primary graph as well as 
implement preparation, manipulation, and read-out. 

Another possibility is to combine a set of recently 
demonstrated experimental capabilities. Wide lattice 
spacings on the order of 5 pm have been achieved [Slj , 
providing a long coherence time to the sites and effec- 
tively approximating the infinitely deep lattice of the 
primary graph. A quantum gas microscope, employing 
a high-numerical-aperture lens, has been used to image 
individual sites in a traditional optical lattice (s^ . Re- 
purposing such a system to focus a laser to a similar res- 
olution would provide a method of manipulating qubits 
by addressing single sites in the case of a Z gate, or mod- 
ifying the potential between sites in the case of X and 
CPHASE gates. Most recently, arbitrary configurations of 
atom positions within a lattice have been implemented 
[ssj , which in particular would allow for the straightfor- 
ward preparation of the initial jO)*^" state as a single 
straight line, one atom wide. In each of these experi- 
ments, read-out of the positions of the trapped atoms is 
also performed. 



V. CONCLUSIONS 

The Bose-Hubbard model, a well-known and widely 
applicable description of bosons confined to a lattice, can 
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be used to generate a universal set of quantum logic gates. 
These gates act on quantum information encoded in the 
positional states of spinless bosons, in contrast to the 
standard method of using internal degrees of freedom 
to store information. Despite arising from Hamiltoni- 
ans that couple outside of the computational space, the 
gates are nevertheless unitary when the evolution param- 
eters are tuned appropriately. The concept of encoding 
a computational space as a subspace of a larger physical 
Hilbert space is often considered to be detrimental or ar- 
tificial when couplings are present between the subspace 
and the larger space. We have shown that such a sys- 
tem can gain power from its ability to access that larger 
space during evolution, while nevertheless remaining re- 
stricted to the computational subspace by the end of its 
evolution. 

Bosons hopping on a lattice under the Bose-Hubbard 
model can also be interpreted as multiple interacting 
quantum walkers on a primary graph that encodes the 
sites and tunneling amplitudes of the lattice in its vertices 
and edges, respectively. The use of multiple quantum 
walkers allows for a sequence of graphs on 0{n) vertices 
to encode a quantum computation on n qubits by dis- 



continuous walk. Such a setup can also be interpreted as 
a single quantum walker on a graph with $7(2") vertices 
and connected to the standard rail model for computa- 
tion by quantum walk, showing the power of multiple 
walkers to eliminate the need for exponentially growing 
resources. Furthermore, localized single-qubit measure- 
ment operators are now possible in a quantum-walk sce- 
nario, and we have presented one example of a scheme to 
implement quantum error-correcting codes with multiple 
quantum walkers under the Bose-Hubbard model. 

Finally, we have discussed the possibility of adapting 
current experimental methods to implement a proof-of- 
principle version of our proposal. Trapped neutral atoms 
in optical lattices with various methods for addressing 
individual sites and manipulating the potential landscape 
fulfill the requirements of such an implementation. 
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